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We present superfiuid density calculations for the unconventional superconductor PuCoGa 5 by 
solving the real-space Bogoliubov-de Gennes equations on a square lattice within the Swiss Cheese 
model in the presence of strong on-site disorder. We find that despite strong electronic inhomogene- 
ity, one can establish a one-to-one correspondence between the local maps of the density of states, 
superconducting order parameter, and superfiuid density. In this model, strong on-site impurity 
scattering punches localized holes into the fabric of d-wave superconductivity similar to a Swiss 
cheese. Already a two-dimensional impurity concentration of nj mp = 4% gives rise to a pronounced 
short-range suppression of the order parameter and a suppression of the superconducting transition 
temperature T c by roughly 20% compared to its pure limit value T c o, whereas the superfiuid density 
p s is reduced drastically by about 70%. This result is consistent with available experimental data 
for aged (400 day-old) and fresh (25 day-old) PuCoGas superconducting samples. In addition, we 
show that the T 2 — dependence of the low-T superfiuid density, a signature of dirty d-wave supercon- 
ductivity, originates from a combined effect in the density of states of 'gap filling' and 'gap closing'. 
Finally, we demonstate that the Uemuera plot of T c vs. p 3 deviates sharply from the conventional 
Abrikosov-Gor'kov theory for radiation- induced defects in PuCoGas, but follows the same trend of 
short-coherence-length high-T c cuprate superconductors. 

PACS numbers: 74.70.Tx, 74.62.En,74.25.N-,76.75.+i 



I. INTRODUCTION 

The study of the superfiuid density or superfuid stiff- 
ness in response to disorder or defects can provide valu- 
able information about the nature of the superconducting 
order parameter. The superfiuid density p s is propor- 
tional to the inverse-square of the penetration depth A, 
a characteristic length of a bulk superconductor deter- 
mining the penetration of a magnetic field inside it. The 
low-temperature dependence of A(T) is being studied ex- 
tensively to unravel the pairing symmetry of the ground 
state of bulk superconductors. In conventional Bardeen- 
Cooper-Schrieffer (BCS) superconductors, with a fully 
gapped excitation spectrum, the penetration depth ex- 
hibits exponential behavior at low temperatures. In 
contrast, in unconventional superconductors, with nodal 
lines or nodal points of the gap function on the Fermi 
surface, A(T) shows power-law behavior depending on the 
type of nodes PQ] Therefore the observation of power laws 
in AA(T) = A(T) — A(0) has been taken synonymous with 
unconventional pairing symmetries in the heavy-fermion 
and high-temperature copper-oxide superconductors. 

It is known that disorder changes the power-law be- 
havior from linear to quadratic in T for d-wave super- 
conductivity as disorder fills in impurity states in the 
nodal gap regions PED Since then it was shown that dis- 
order in unconventional superconductors leads to an even 
stronger suppression of the superfiuid density p s rela- 
tive to its superconducting transition temperature T c 
than expected from unitarity scattering in the Abrikosov- 
Gor'kov (AG) theory. In particular, Franz and co- 
workers^ investigated short-coherence-length supercon- 
ductivity in cuprates, that is, the coherence length £ is 
comparable to the lattice parameter a, and demonstrated 



that it behaves markedly differently from the AG the- 
ory of impurity-averaged order parameters. In fact, de- 
tailed calculations of the spatial dependence of the lo- 
cal density of states and order parameter in the vicin- 
ity of an impurity with strong (unitarity) nonmagnetic 
scattering potential showed that the order parameter is 
abruptly suppressed within just a few lattice parame- 
ters resembling the holes within a Swiss cheese EaED This 
is in stark contrast to the assumption in the AG the- 
ory of pair-breaking with a spatially uniform suppressed 
impurity-averaged order parameter. The Swiss Cheese 
model was originally introduced to explain the univer- 
sal scaling behavior of superconducting transition tem- 
perature versus zero-temperature superfiuid density in 
the Uemura plot, T c vs. p s (0)p ^ of underdoped high- 
Tc cuprate superconductors^^]] anc j found a recent re- 
vival for describing the behavior of Kondo holes in heavy- 
fermion superconductors.^ 

In this work, we use the Bogoliubov-de Gennes 
(BdG) lattice model to study the effects of disorder on 
bulk and local properties in superconducting PuCoGas. 
Plutonium-based superconductors can be thought of as 
superconducting clocks, since the natural radioactivity 
of Pu ( 239 Pu half-life = 24, 000 years) creates lattice de- 
fects like clockwork, which scatter electrons and break 
Cooper pairs. Pair-breaking creates impurity bands in 
the nodal regions affecting the signature of a pure d- 
wave superconductor. Such pair-breaking effects are evi- 
dent in the suppression of T c with timeJ^HUl the spin- 
lattice relaxation rate 1/TiJ^ where the T 3 tempera- 
ture dependence switches to a linear-in-T behavior, and 
the penetration depth A(T), where the linear-in-T tem- 
perature dep enden ce gives way to a T 2 behavior at low 
temperatures P3HD Unlike other bulk probes, the magni- 
tude of A(0) is very sensitive to defects owing to its nature 
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of measuring the stiffness of the superconducting phase 
coherence in the sample, whereas the magnitude of T c is 
less sensitive to defects, since it is related to the spatial 
average of the order parameter. A detailed account of 
this difference will be presented. 

The main result of our study is that self-irradiation- 
induced defects in PuCoGas violate the AG theory of 
dilute disorder in superconductors, which is based on the 
premise of an impurity-averaged order parameter, while 
it is fully consistent with the Swiss Cheese model of short- 
coherence-length superconductors. Our detailed calcula- 
tions show that measurements of the transition temper- 
ature and superfluid density in fresh (25 day-old) and 
aged (400 day-old) PuCoGas are fully consistent with 
the Swiss Cheese model of disorder, which captures both 
the 20% suppression of the superconducting transition 
temperature T c and a large 70% suppression of the su- 
perfluid density p s (0) between aged and fresh samples, 
as well as the temperature dependence of p s {T). Fi- 
nally, our study of the Swiss Cheese model exemplifies the 
common behavior of the superfluid density in unconven- 
tional short-coherence-length superconductors with dis- 
order. Because of spatially extended quasiparticle excita- 
tions along the nodal directions of the order parameter, 
in the presence of a strong on-site impurity potential, 
one finds the suppression of p s (T) to extend over several 
coherence lengths, whereas the suppression of the order 
parameter is very localized and limited to a few lattice 
sites. In that respect a close relationship between the 
copper-oxide and plutonium based superconductors ex- 
ists. 

The paper is organized as follows: In Sec. II, we in- 
troduce the BdG lattice formalism of the Swiss Cheese 
model and the corresponding superfluid density expres- 
sion. The local variation of the density of states, super- 
conducting order parameter, and superfluid density for 
a single impurity, as well as for multiple impurities, are 
presented in Sec. Ill A. The evolution of the spatially 
averaged superfluid density and corresponding density of 
states as a function of impurity concentration is discussed 
in Sec. Ill B. In Sec. Ill C we present the Uemura plot 
of T c vs. p s (0) for the Pu-based compound and compare 
with the results for cuprate superconductors. Finally, we 
conclude in Sec. V. 



A. The BdG lattice model 

It has been shown that superconductivity in PuCoGas 
is unconventional because thermodynamic, transport, 
nuclear magnetic resonance, and neutron scattering data 
are consistent with a model based on e?-wave pair- 
ing symmetr y (that is, based on nodal lines in the 
gap function) jmmUHH] First-principles calculations^ 2 ^ 
have shown that PuCoGas hosts a Fermi surface with 
four sheets, of which two are cylindrical sheets centered 
at the M point in the Brillouin zone. The quasi-two- 
dimensional nature is related to the layered structure of 
Pu atoms forming a square lattice. For simplicity, we 
consider only a single-band model, which captures the 
essential physics of <i-wave superconductors. In order 
to describe the disorder effect on superconductivity, we 
consider a tight-binding model Hamiltonian defined on a 
square lattice: 

ijcr icy 



Here c| CT (cjo-) creates (annihilates) an electron at the i-th 
site of spin a. The variables and fi are the hopping 
integrals and chemical potential, respectively. We model 
the disorder by considering a distribution of short-ranged 
nonmagnetic impurities, that is, ej = Ui mp 6u, with Ui mp 
representing the potential scattering strength. The quan- 
tity Ay denotes the superconducting order parameter or 
gap function. Since the origin of superconductivity is not 
our concern, we introduce an effective nearest-neighbor 
pairing interaction V , such that the superconducting or- 
der parameter is determined self-consistently: 



Ay = — (c it c u - c 4 c it ) , 



(2) 



where (ij) is a nearest-neighbor (NN) site pair, and zero 
otherwise. By using the Bogoliubov transformation 



H In ~ V t 7T] , 
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II. THEORETICAL MODEL 

We begin with the BdG mean-field theory of the attrac- 
tive Hubbard model which was used extensively to de- 
scribe superconductivity in correlated electron systems. 
Our main interest is focused on the linear response cal- 
culation of the superfluid density p s for the BdG lattice 
model following the approach described by Scalapino et 
alW^ for the Hubbard lattice model. This generic ap- 
proa ch was consequen tly ap plied to two-dimensional d- 
wav d 12 l 3( ^J and s-wave^2HS! superconductors. 



the Hamiltonian in Eq. (TH) can be diago nalized by solving 
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the corresponding BdG equati 
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Here (m",w™) t are the eigenfunctions at site i corre- 
sponding to the quasiparticle excitation energy E n , and 
the normal-state single-particle lattice Hamiltonian is 



H 



(5) 
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Throughout this work, we limit the hopping integrals 
only to the NN sites on the square lattice, that is, 
U,i+s = t for 6 = (±1, 0) and (0, ±1), and zero otherwise. 
The self-consistency equation for the superconducting or- 
der parameter on the square lattice is thus given by 



V 
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(6) 



where S = (±1,0) and (0,±1), and zero otherwise, and 
the temperature is denoted by T. For d x 2_ y 2-wave pair- 
ing symmetry, the order parameter along the y direction 
has opposite sign compared to the x direction, which is 
indeed obtained in the solutions. Once the BdG equa- 
tions are solved, many interesting properties can be ex- 
plored. For example, for a given concentration of disor- 
der, the superconducting transition temperature can be 
determined by the condition that the averaged supercon- 
ducting order parameter vanishes. Another important 
observable that can be tested experimentally is the local 
density of states (LDOS) at zero temperature, which is 
given by 

Pl (E) = 2Y / lK\ 2 6(E-E n ) + K\ 2 S(E + E n )] , (7) 



K x {ri) = - ^2[ti. l+ scl a c i+ s, a + h.c] , (11) 

are the particle current and kinetic energy densities. The 
variable S — x,x±y denotes the links which have contri- 
bution to the bond current and kinetic energy along the 
x axis. The charge current density operator along the x 
axis is then found to be 
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SH'(t) 
SA x (r h t) 



eJf (r<) + e 2 K x {vi)A x {v u t) . 

(12) 

We calculate the paramagnetic component of the electric 
current density to first order in A x , 



(Jf (r s: )) = -i / {[jF{t),H'{t')]-) dt', (13) 



and the diamagnetic part in {K x )q only to zeroth order; 
(. . . )o represents a thermodynamic average with respect 
to Hq. Straightforward algebra yields the current re- 
sponse function 



where the prefactor '2' is due to the twofold spin degen- 
eracy. The differential tunneling conductance in scan- 
ning tunneling spectroscopy is directly proportional to 
the LDOS and can provide insights into the electronic 
properties and symmetry of a superconductor. Since one 
of the major pieces of interest in the present work is the 
superfluid density, which characterizes the superconduct- 
ing phase rigidity, it is useful to give an expanded dis- 
cussion of its derivation and lattice formulation. Our 
BdG lattice formulation of the superfluid density follows 
closely the seminal work by Scalapino and coworkers^ for 
the Hubbard model on a lattice. We calculate the super- 
fluid stiffness for a current response to a vector potential 
of wave vector q and frequency uj along the x direction as 
given by the Kubo formula. For this purpose, we expand 
the Hamiltonian to include the interactions of electrons 
coupled to an electromagnetic field. The tile-dependent 
total Hamiltonian is 

H t = H a + H'(t). (8) 

Here H'(t) describes such a minimal coupling 

H'(t) = -ea J^M^t) (^N + y^( r i>*)^( r i)) - 

i 

(9) 

where a is the lattice constant, A x is the vector potential 
along the a; axis, and 

Jx( r i) = -i^iUA+s^Ct+s^ - h - c -] > ( 10 ) 

a.8 



e 2 A x (r. l ) 



dt'([jf(q,i),4 P (-q,0]-)o 

1 

■{KxMh . (14) 



By performing a lattice average over the variable r.; to 
eliminate the atomic-scale fluctuations, we define an ef- 
fective "Drude weight" as a measure of the superfluid 
density 



Ps= ^ = -<iQ+n KX (q^0,u; = 0) . (15) 



The first term is the kinetic energy along the x direction 
divided by the number of lattice sites, 



+(l-f(E n )M +Xia vf:+c.c]\ . (16) 



It represents the diamagnetic response to an external 
magnetic field B — V x A with gauge A. x ^ and 
A y = A z = 0. The second term is the paramagnetic re- 
sponse given by the disorder-averaged transverse current- 
current correlation function 
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with coefficients 

A ni ,„ 2 (q) = 2^e-^[<5< 2 -<^<^],(18) 

i 

D ni , n2 (q) = 2^e-^[<>r-«r<5]-(19) 

i 

B. Numerical solution of the BdG equations 

We follow an iterative numerical procedure to solve 
self-consistcntly the BdG equations via exact diagonal- 
ization: For a given impurity distribution, we start with 
a uniform d 2 .2_ y 2-wave order parameter, that is A^j+i = 
— Aj ; j_££, at a low temperature. After the BdG equa- 
tions Q are diagonalized, the obtained eigenvalues and 
eigenfunctions are used to update the superconducting 
order parameter as given by Eq. Then we start a 
new cycle of iteration. The iteration will continue until a 
convergence criterion for the order parameter is reached. 
We have taken the difference of the order parameter at all 
sites between two consecutive iterations to be less than 
10~ 8 as the convergence criterion. Upon convergence is 
reached at a given temperature, both the local density of 
states (LDOS) and superfluid density are calculated by 
using the supercell technique to reduce finite-size effects. 
The final LDOS and superfluid density are the result of 
an ensemble average over about 20 impurity configura- 
tions. After this step, the calculation is moved to the 
next higher temperature point. 

Next we also comment on the choice of lattice model 
parameters. In our numerical calculations, we have cho- 
sen to measure the energy in units of the nearest-neighbor 
hopping parameter t with chemical potential /i = — 0.36i 
and superconducting pairing interaction strength V — 
1.13i. The on-site impurity potential scattering strength 
Ui mp = lOOt was used to model the strong (unitarity) 
limit of impurity scattering. Additionally, we broadened 
the density of states calculations by a small imaginary 
term of width T — O.Oli to overcome the discreteness of 
the energy spectrum. The size of a single cell for the 
self-consistent lattice model calculation is at least of 20- 
by-20 sites, when the supercell method with 6-by-6 cells 
is employed or 35-by-35 sites otherwise. However, we still 
encountered finite size effects for systems of such size at 
low temperatures and for low disorder concentrations. 



III. RESULTS AND DISCUSSION 

A. Real-space imaging of local properties around a 
single impurity 

We begin with spatial images of various local proper- 
ties at T = 0, including the LDOS at zero excitation 
energy Nq(t), the superconducting (SC) gap amplitude 
|A(r)|, and the superfluid density p s (r) all shown in the 
top row of Fig. 1 for a single impurity at the center of 
the cell. The SC gap amplitude at a particular site (x, y) 
is obtained from | A(r)| = [A(x + Sx, y) + A(x — Sx, y) — 
A(x,y + 5y) — A(x,y — <%)]/4, where S is the distance 
between NN sites. [Below we denote A(0) as the average 
SC gap at T = 0.] Whereas p s ( r ) is the current-current 
correlation function which involves a double summation 
over real-space. Inserting Eqs. (18)-(19) into Eq. 18, we 
obtain p s (r) ~ J2i j ^(A^ with site indices 

We plot p s (r) in Fig. l(cl) to visualize the effect of a 
single impurity on the rigidity of the superconducting 
phase. 

The well-known Friedel oscillation is clearly seen in 
the LDOS Nq(t) along the nodal direction of the SC 
gap. This is a signature of unconventional supercon- 
ductors, whi ch wa s studied extensively in the high- 
T c cuprateJ ^ 16 * 18 * and iron-based superconductors!^^ 
The d x 2_ y 2-wave order parameter is fully suppressed at 
the impurity site for a strong scatterer, aside from very 
weak oscillations near the order parameter maximum, see 
Fig. l(bl). The spatial dependence of the superfluid den- 
sity p s (r) is also strongly localized around the impurity 
with characteristic features more similar to the LDOS 
than the order parameter. To illustrate this type of Swiss 
cheese phenomenon of strongly suppressed superconduct- 
ing properties near an impurity, we have taken various 
cuts along (100), (010), and (110) directions, as shown 
by solid red, dashed blue and solid green lines, respec- 
tively, in the middle row in Fig. 1. 

Next we revisit characteristic properties of a single 
impurity in a <i-wave superconductor. The LDOS is a 
measure of the single-quasipart icle s pectral weight. In 
agreement with earlier works^ H 16 * 18 ! it decays inversely 
quadratic with distance r from the impurity, 1/r 2 , along 
the nodal direction of the gap (green line), whereas 
along the antinodal direction, the decay is exponential. 
The spatial order parameter amplitude |A(r)|, plotted 
in Fig. I(b2), creates a resonance state at the impu- 
rity site, reflecting the pair-breaking characteristics of 
d-wave pairing. The fourfold modulation of the order 
parameter is preserved in the case of a scalar impurity, 
while the suppression is smoother along the nodal di- 
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FIG. 1: (Color online) Visualization of the local behavior of superconducting properties for a single impurity at the center of the 
cell, see panels (al)-(cl). The 3D plots of No(r) in (al), A(r) in (bl) and p s (r) in (cl) show the real-space modulation of these 
properties in response to a single impurity at the center. The corresponding color-coded contour maps highlight the patterns of 
the modulations. All spectra are calculated on a 35-by-35 lattice and then interpolated for visualization. The inset to (al) gives 
the schematic view of various directions and locations of sites with respect to an impurity (red dot) at which the following plots 
are drawn. Each plot in the middle row corresponds to a ID cut through the spectrum shown in the corresponding top panel 
(the color of each representative curve is the same as for the arrow drawn in the inset). The shading at the lattice boundaries 
delineates the region where finite-size effects are expected to affect the results shown in panels (a3)-(c3). Similarly the curves 
in the bottom row are drawn at four representative sites with respect to the location of the impurity and are compared with 
their average value. Site 'A' is the nearest neighbor to the impurity along the antinodal direction. The next-nearest-neighbor 
site 'B' is along the nodal direction. Sites 'C and 'D' are farther away from the impurity location along the antinodal and 
nodal directions, respectively. 



rection where strongly damped Friedel oscillations are 
found. In contrast, the superfiuid density p s ( r ) shows 
the usual fourfold modulation but in addition picks up 
the phase of the d x 2_ y 2-wsve pairing symmetry as shown 
in Fig. I(c2). Along the (lOO)-antinodal direction, where 
the gap is positive, p s shows a remarkable enhancement 
from its average value at the NN lattice site from the 
impurity (red line), whereas the NN lattice site along 
the (OlO)-antinodal direction exhibits stronger suppres- 
sion (blue dashed line). This result is expected in linear 
response for a supercurrent flowing along the (100) di- 
rection, because it breaks the tetragonal symmetry of 
the lattice. Along the diagonal direction, the power-law 
behavior of p s arises from gapless quasiparticles in the 



LDOS JV (r)PE5l 

To gain further insight into the impurity effect on p s (r) , 
we plot the temperature evolution of these quantities at 
four representative sites with respect to the impurity po- 
sition and compare them with the average value. In 
the case of a clean superconductor with lines of nodes 
in the gap function, the low-temperature (fc^T <C A) 
approximation yields a linear behavior of the superfiuid 
density for gap less nodal quasiparticles on the Fermi 

surfaC ePn™i 

p s (T)/p s (0) ~ 1 - dT. (20) 

Here, the slope C\ is determined by the BCS ratio 
A(0) /T c . In the presence of an impurity with large broad- 
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ening ^> fc^T, the lowest order temperature dependence 
of the superfluid density becomes T 2 [Refs. I7HTT] (note 
that for an isotropic s-wave gap the low-energy superfluid 
density decays exponentially): 

p s (T)/ Ps (0)~l-C 2 T 2 . (21) 

Here C'2 is a more complex function of Fermi velocity 
and effective impurity broadening.^ I n Fig. I(c3) we see 
that at the NN site 'A' along the (100) direction, where 
p s (A) shows a peak, the corresponding p B (A,T) shows a 
quasilinear behavior at low temperatures, within the er- 
ror produced by finite-size effect. At the NNN site 'B' 
along the node, the superfluid density is suppressed show- 
ing a T 2 dependence in accord with the average value. 
We believe that in principle magnetic field-angle depen- 
dent measurements of the magnetization can provide a 
unique and indispensable tool to probe the presence of 
gap nodes and shed light on the pairing symmetry in the 
ground state of a bulk superconductor. This idea is re- 
lated to the observation o f the nonlinear Meissner effect 
predicted by Yip and Sauls^J^I f or fields along the nodal 
and antinodal directions. However, its analysis will be 
complicated by surface bound states.^ 

Next we turn to the case of randomly distributed im- 
purities in Fig. 2 for a concentration rii mv = 2%, where 
the local information of the previously discussed quan- 
tities changes more dramatically across the entire lat- 
tice. Nevertheless, the local response of all quantities 
with respect to the impurity location as well as the one- 
to-one correspondence between LDOS, order parameter, 
and superfluid density is present for all impurity concen- 
trations studied. Differences between multiple impurities 
and single impurity are appreciable in all calculated spec- 
tra. Due to the quasiparticle interference and the overlap 
of wavefunctions of quasiparticles scattering off impuri- 
ties, Friedel oscillations are enhanced and present for all 
directions shown in Fig. 2(a2)-(c2). Similarly, the dif- 
ferences in the T dependence of properties at each site 
is more clearly visible here as can be seen in Fig. 2(a3)- 
(c3). For example, at position 'B', which sits at the cen- 
ter along the nodal direction between two impurities, one 
can probe the gapless quasiparticles and thus show quasi- 
linear behavior in p s (T) at low T. Whereas at site 'A', 
which sits nearly at the center between two impurities, 
but along an off-nodal direction, one sees enhancement 
in A(0) and p s at low T. 



B. Temperature dependence of average superfluid 
density 

Figure 3 presents the calculated DOS (averaged over 
the entire 35-by-35 lattice sites) as a function of energy 
for the d-wave pairing case and impurity concentrations 
nimp — 0, 1, 2, 3, 4%. The disorder-induced change in the 
DOS is also reflected in the temperature dependence of 
the spatially averaged superfluid density. Since finite-size 



effects are more pronounced at low energy and for low im- 
purity concentrations, for example, the zero-energy DOS 
of the pure sample (red line) in Fig. 3(a) is expected to 
vanish, while it is of the order of the numerical broaden- 
ing term ~ >/r/A(0) N n (0), where N n (0) is the normal 
state DOS at the Fermi level. For the same reasons, the 
low-temperature values of p s (T) are expected to be less 
accurate. Nevertheless, our calculations reproduce the 
hallmark W-shaped feature of the DOS due to d— wave 
pairing for all concentrations. It is interesting to point 
out that, with increasing concentration not only is the 
gap amplitude decreasing, but also the gap nodes be- 
come more filled in by impurity states, see shaded and 
colored areas in Fig. 3(a). The gap-filling or build-up 
of resonant impurity states near the Fermi level (which 
can be quantified by an impurity scattering rate) are 
the main aspect of impurity effects that give rise to the 
quadratic-in-T dependence of the superfluid density at 
low temperatures. For large impurity concentrations the 
Swiss Cheese model recovers the dirty d-wave result of 
Eq. (21) for the temperature dependence of p s (T). Plot- 
ting p s (T) as a function of T 2 in the inset of Fig. 3(b) 
demonstrates the gradual change from linear to quadratic 
behavior with increasing impurity concentration. 

To connect our Swiss Cheese model calculations with 
muon-spin rotation measurements of the penetration 
depth, we compare in detail in Fig. 3(c) the temperature 
dependence of p s {T) of the fresh and aged samples with 
our BdG calculations in the limit of strong on-site impu- 
rity scattering. The impurity concentrations ni mp = 1% 
and 4% were chosen to reproduce the observed suppres- 
sion of T c and p s (0) of the fresh and aged samples, re- 
spectively. The relevant model parameters were given 
in Sec. |II B| which are related to the coherence length 
£ ~ hvp/ (7rA), where the Fermi velocity is given approx- 
imately by v F ~ at/h with A = 4A(0) and A(0) = 0.08i 
at n imp — 4%, so that £ ~ 3.9a and lattice parame- 
ter a = 0.423 nmP From our calculated T c suppression 
and comparison with experiment, we conclude that hy- 
pothetically pure PuCoGas has a bare superconducting 
transition of T c q — 18.9 K in agreement with previous 
estimates. The reported decrease in T c of about 3 K 
is in agreement with the radiation-induced reduction of 
T c (ps 0.24 K/month) reported for PuCoGas samples of 
slightly different isotopic concentrations J5HE5I We list in 
Table I all characteristic parameters. Measurements on 
the fresh sample (after 25 days) were performed with an 
applied field H = 60 mT. The aged sample (after 400 
days) was measured in applied fields of 60 mT and 300 
mT. All measurements were performed in a field-cooled 
mode above the lower critical field H C \W^ 

The excellent agreement between the short-coherence- 
lcngth BdG calculations within the Swiss Cheese model 
and the measured T c and superfluid density p s {T), com- 
bined with the previously demonstrated failure of the 
dirty d-wave theory, 27 shows that the uniform, dilute- 
impurity pair-breaking theory by Abrikosov and Gor'kov 
is not applicable to PuCoGas. 
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FIG. 2: (Color online) Visualization of the local behavior of superconducting properties for a randomly distributed impurity 
concentration of 2%, same notation as in Fig. 1. Brown cross symbols on the 2D color maps in the top row depict the locations 
of impurities on the lattice. The selection of test sites used in the bottom row of panels is the same for all spectra and marked 
in panel (al). Unlike for the single-impurity case in Fig. 1, the SC gap and superfiuid density now show much stronger spatial 
fluctuations. All averaged calculations are done for 10 samplings over random configurations of impurities. The visualized 
results are for one particular disorder configuration. 



TABLE I: Parameters derived from comparing A(T) in 25 day-old (fresh) and 400 day-old (aged) PuCoGas to the BdG 
lattice model within the Swiss Cheese model of d^.^-wave superconductivity, penetration depth A(0) (at T = 0), impurity 
concentration ni mp , and the transition temperature of the nominally pure sample T c q ~ 18.9 K. 



Samples 




T c in K 


T c in K 


A(0) (nm) 


A(0) (nm) 




[Present Theory] 


[Present calculation] 


[From Refs. 125126] 


[Present calculation] 


[From Refs. 125(26] 


Fresh (25 days) 


1% 


18.88 


18.25(10) 


310 


265(5) 


Aged (400 days) 


4% 


16.45 


15.0(1) 


524 


498(10) 



C. Suppression of p s (0) and T c 

The magnitudes of the superconducting properties 
N(0), T c and p s (0) provide valuable information about 
the topology of the gap function on the Fermi surface as 
a function of impurity concentration, see Fig. 4(a). Fur- 



thermore, the correlation between them exemplifies the 
deviations from conventional dirty d-wave theory arising 
in a system with short-coherence-length superconductiv- 
ity. While the dependence of N(0) and A(0) (also T c ) on 
disorder is quasilinear, the suppression of the superfiuid 
density is dramatically enhanced. These features demon- 




FIG. 3: (Color online) The energy dependence of the DOS is related to the temperature dependence of p s (T) (averaged over 35- 
by-35 lattice) for different impurity concentrations, (a) Computed DOS N(u>) is plotted for increasing impurity concentration 
along the vertical axis (curves are not shifted vertically). The grey shaded area highlights the nature of the gap closing, while 
the colored filling in each spectrum illustrates the trend of gap filling with increasing disorder. Consequences of finite-size effects 
can be seen in the value of the zero-energy DOS for the pure sample (red line), which is expected to be zero in the absence of 
a numerical broadening term F. (b) Calculated results of p s {T) (normalized to their corresponding zero temperature value). 
The inset expands the low-T region of p B (T) vs. T 2 to emphasize the quadratic temperature dependence, (c) Computed p B (T) 
for 1% and 4% impurity concentrations are compared with the measured data for a fresh (25 day-old) and old (400 day-old) 
PuCoGas sample, respectively.^ The cyan diamond and red circle symbols are for field-cooled measurements in Ho = 60 mT, 
while the cyan star symbols are for Ho = 300 mT. Since no data are available at zero temperature and near T c , we extrapolated 
each experimental curve into these regions for proper normalization. After extrapolating the data, we estimated the values of 
p s (0) and T c , which are slightly higher than the ones predicted using a linear extrapolation methodp^' see Table I. 



strate the strong deviation of the Swiss Cheese model 
from conventional AG theory as discussed in details be- 
low. It has been noted before that the suppression of the 
zero-temperature superfiuid density for d-wave pairing 
with disorder on a square lattice is much stronger th an for 
the superconducting gap or transition temperature! 12 * 30 ^ 
In Fig. 4(b) we show that A(0) and T c exhibit a weakly 
linear decrease with an increasing DOS at the Fermi level, 
N(0). On the other hand, p s (0) decreases faster than 
linear with increasing N(0). The superfiuid phase coher- 
ence is destroyed more rapidly than the superconducting 
amplitude accounting for the marked difference between 
fresh and aged samples. The dependence of the super- 
fluid density on the order parameter is given in Fig. 4(c), 
while the same information as a function of T c is pre- 
sented in the Uemura plot in Fig. 5. The BCS ratio 
A(0)/T c increases gradually with increasing gap ampli- 
tude, that is, with decreasing impurity concentration as 
expected. As mentioned before, p s is suppressed much 
faster than T c or A(0) For ni mp = 4% we calculate for 
T c or A(0) a suppression of ~ 20%, whereas p s (0) is sup- 
pressed by almost 70%. This finding is consistent with 
the experimental data for cuprates where the supercon- 
ducting transition temperature is much more robust to 



disorder than what is predicted by the AG theory for d- 
wave pairing, when measured against the corresponding 
change in p s (0). 

To elaborate some more on this point, we plot the ex- 
perimental data of p s (0) vs. T c in Fig. 5 for PuCoGas 
by Ohishi et a/p3 and compare with available experi- 
mental data on YBa 2 Cu307_,5 (YBCO) superconductors, 
the self-consistent T-matrix result of the dirty d-wave 
theory,^ and the Swiss Cheese model of the BdG lattice 
model. The situation of PuCoGas is not unprecedented. 
For example, in YBCO p s (0) is suppressed much more 
dramatically than T c . In Ni-doped and He-irradiated 
YBCO, a suppression of T c by about 20% is accompa- 
nied by a suppression of p a (0) by about 70% PSEEM 
Several years ago, the same underlying disorder physics 
has been discussed within the BdG lattice modelP ^ 60 * 61 ! 
Where applicable we find agreement with these calcu- 
lations. Franz et aP^I compared their BdG results to 
YBCO samples for various conditions of 

disordering 

and concluded that the effect is enhanced by a short co- 
herence length £/a «2-5p^ As in YBCO, the coherence 
length in PuCoGas is relatively small (~ 2 nm in both 
materials) with comparable lattice parameters (~ 0.5 nm 
in both materials). Note that a similar enhanced sup- 
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FIG. 4: (Color online) Disorder induced correlations between 
N(0), A(0), T c and p s (0). (a) Variation of different properties 
(see legend) as a function of impurity concentration rii mp %. 
(b)-(c) The same quantities are drawn as a function of total 
number of quasiparticles on the Fermi surface in (b) and as a 
function of order parameter suppression. 
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FIG. 5: (Color online) Uemura plot of superfluid density in a 
disordered system. The two red open squares are the experi- 
mental data for PuCoGa 5 (Pu-fI5) from Ref. [27] The rest of 
the open symbols are taken from YBCO data under different 
conditions of disorder environments!^ Data for YBCO films 
are obtained from Ref. I54I for YBCO crystals from Ref. I55I 
for ceramic samples from Ref.[56]and for He-irradiated YBCO 
films from Ref. I57I The filled red dots are the present BdG 
theory. 



pression of the superfluid density with chemical doping 
has been reported for the related cerium-based compound 
CeCoIns-zSna;, where £ ~ 3 nmP 

Taking into account the considerable spread among all 
data sets, a fit gives T c /T c q oc [p s (ni mp ) / p s {0)] 0A , which 
deviates significantly from the linear scaling of the Ue- 
mura plot (T c oc p s (0)) of underdoped high-T c cuprates, 



which has its origin in strongly correlated electron inter- 
actions. Finally, it is evident that the AG theory over- 
estimates the suppression of T c , while the Swiss Cheese 
model is in very good agreement with all data sets. 



IV. CONCLUSIONS 

We have found good agreement between the results of 
the Swiss Cheese model using the BdG lattice model and 
the superconducting properties of PuCoGas. Most im- 
portantly, the results demonstrate that, despite strong 
electronic inhomogeneity, a one-to-one correspondence 
between the electronic and superconducting linear re- 
sponse functions on each lattice site is maintained at all 
impurity concentration. 

The low-temperature dependence of the superfluid 
density or penetration depth in both fresh (25 day-old) 
and aged (400 day-old) PuCoGas are consistent with a 
line of nodes in a strongly disordered d x 2_ y 2--wave order 
parameter. The Swiss Cheese model can describe the 
quadratic temperature dependence of the superfluid den- 
sity p s (T) with the gap A(0) or T c reduced by about 20% 
for impurity concentration nj mp = 4%. It can account for 
at least a 70% reduction in p s (0) contrary to the dirty d- 
wave theory. We attribute this to the fact that PuCoGa 5 
possesses a relatively short coherence length, and, there- 
fore, the conventional Abrikosov-Gor'kov pair-breaking 
theory, in which the order parameter is spatially aver- 
aged, is inappropriate. This result is similar to what is 
known for radiation damaged or doped high-T c cuprate 
superconductors. 

For low impurity concentrations the short-cohercncc- 
length model agrees with the expected suppression of 
p s {0) obtained from the dirty d-wave theory, while it de- 
viates drastically otherwise. Irrespective of the approach 
to impurity averaging, both theories predict a change 
from linear to quadratic in temperature in the superfluid 
density for large impurity concentrations. We show that 
the quadratic dependence of the superfluid density arises 
from the combined effects of the SC 'gap filling' and 'gap 
closing' in response to the presence of strong disorder. 

Furthermore, both the fresh and aged PuCoGas sam- 
ples are consistent with the Swiss Cheese model for a 
weak-coupling gap A/T c ~ 1.6 - 2.0 [A = 4A(0), aver- 
aged over four NN sites for d-wave pairing] that is sup- 
pressed by strong impurity scattering. Finally, our lattice 
calculations show that although the order parameter is 
significantly suppressed in the immediate vicinity of im- 
purities, and the superfluid density is strongly suppressed 
over extended regions along the nodal directions, the su- 
perconductivity remains remarkably resilient. These cal- 
culations provide further evidence that PuCoGas is the 
link between low-temperature heavy-fermion and high- 
temperature cuprate superconductors. 
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